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Tractor motors always operate in the speed region higher than rated speed, 
but is limited to the module of the stator current, stator voltage vectors. 
Additionally, mathematical model of traction motor has shown nonlinearity 
through the product of the state variables isa, isọ with the input variable 
Ws:Wsisq, Wsisq. Therefore, this paper focuses on the study of speed control 
of traction motors in weakening field region while optimizing torque control, 
and choosing the backstepping method in designing speed—flux controller in 


order to solve the nonlinear structure. The simulation results of the 
responses: speed, torque, power, and flux performed on MATLAB/Simulink 
software with parameters collected from metro Nhon-Hanoi Station, 
Vietnam have proven the correctness in theoretical research. 
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1. INTRODUCTION 

With the outstanding advantages of the ability to carry a large number of passengers, safely, on 
time, reducing traffic congestion and pollution environment in big cities, urban electrified trains are 
increasingly attracted the attention from researchers and governments all over the world [1], [2]. In Vietnam, 
some urban railway lines have been put into operation and others are about to be operated in the near future. 
However, to ensure the quality of train operation, studies on traction engines are essential. Due to the great 
development of power electronic converters, induction motors are widely used as traction motors due to their 
simple structure and speed control range. The operation characteristics of traction motors used in 
transportation different from industrial loads, which is often operated in the speed region greater than the 
rated speed. However, when the control speed is higher than the rated speed, there are limitations on stator 
voltage and current modules, so it is imperative to control spedd in weakening flux region. 

Mengoni et al. [3], Harnefors et al. [4] focused on the stator flux vector control of induction motor 
(IM) in the field weakening region to ensure maximum torque but ignore the requirements for maximum 
torque output. Son et al. [5] proposed a flux control method for operating IM with the condition that voltage 
on DC-link varies to ensure maximum output power. Zarri et al. [6], Abu-Rub and Koltz [7], Shin et al. [8], 
Casadei et al. [9] also studied IM control in weakening flux region with optimal flux selection and control 
methods. IM control by stator flux oriented (SFO) method with stator flux inversely proportional to motor 
speed. However, field oriented control method (FOC) with optimal flux levels clearly calculated as a function 
of machine parameters; controlled by direct torque control (DTC) method with controllers being PI. Tang and 
Huang [10] designed the speed control of IM in weakening field region by regulating PID parameters. 
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Gallegos et al. [11] proposed a new current controller making it easier for IM to operates smoothly and 
maximum torque-per-ampere (MPTA) control in weakening field region given variation of motor parameters 
and voltage source. Casadei et al. [12] presented three sensorless control schemes for IM operation in the 
weakening flux region with maximum torque utilizing maximum voltage and current. 

Lin et al. [13] proposed the backstepping adaptive motion control method to determine the slip 
surface for the sliding controller applied in motion tracking of the mechanical transmission system. Yu et al. 
[14] control position tracking with uncertainty and load disturbance using adaptive fuzzy backstepping 
controller. Echeikh et al. [15] designed non-linear backstepping controller for five-phase induction motor 
drive system operating at low speed. Regaya et al. [16] have designed 2 steps: first step of evaluating rotor 
flux based on the backstepping observer and the second step determines the adaptive mechanism of the rotor 
resistance based on the estimated rotor flux. Kirad et al. [17] applied sensorless backstepping technique with 
kalman filter to improve performances of motor. Belkacem et al. [18] made comparison between 
backstepping and backstepping sliding mode controller for some components of vehicles to show efficiency 
and robustness against external disturbances and various road profile. Nguyen et al. [19], Do et al. [20] used 
backstepping method for IM with constraints of input and output. 

However, the above works have not proposed a combination of nonlinear control methods for 
traction motor-IM operating in the weakening field region applied to means of transportation. Therefore, this 
paper presents speed control of traction motor for urban electrified train in field weakening based on 
backstepping technique applied to metro Nhon—Hanoi Station, Vietnam. 


2. MODELING ELECTRIC TRAIN 

In modeling electrified train, modeling traction motor and load very needs. The continuous state 
model of traction motor is designed on dq coordinate system, analyzing the forces that impede motion 
including friction, air resistance, slope resistance, and curve resistance is shown in Figure 1. Survey of 
parameters and forces conducted on the urban railway line from Nhon to Hanoi Station, Vietnam. 
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Figure 1. Forces act on train 


2.1. Modeling traction motor 
The traction electric motor is induction motor with the field oriented control method based on the dq 
coordinate system. The mathematical are expressed as (1) [21], [22]: 


ora = A + Jisa + Wsisg + oa ra + owa $ wig sa 

ae = Ga t oia T Ossa + Ge Wa =g Ora + i, sn 

Hera isa pra + (ws — w)Prq (1) 
a pisa mm a 7 ra 

My = iZ Dt aksg = ai — bsPratsg 

My = mM, + 2 
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; n Ln Ls . 
Where parameters are defined: Y,a = Wra/Lm Va Vig = Wra/Lm; 6 = 1- re is leakage factor; T, = - is 
; Lps . oL Lin i 
stator time constant; T, = — is rotor time constant; T} = —2,@, = w + Z% sa. 
Rr To Tr Wra 
is pairs of poles; J is torque of inertia of electric train and motor; m, is load torque; Yra, Prq is rotor flux; 


Lm Lr» Ls are mutual, rotor, and stator inductances respectively. Inertia torque of electric train calculated [21]: 


w is mechanical rotor speed; Zp 


Jeq = na (Py Q) 


CAN T 
where Mmas is mass of train, N is numbers of motors, Dwn is diameter of wheel, and t is transmission ratio. 


2.2. Modeling load of train 
The forces impeding the movement of the train include the air resistance, the resistance due to 
rolling friction, which is reduced to the basic resistance by mass calculated based on the Davis formula [23]-[25]: 
wo = — =a +bv + cv? (3) 


Mmas 


Where a, b, c are the drag coefficients provided by the manufacturer and v is the speed of the train. 
Resistance due to slope: 

Foraa = mg sina (4) 
Where sin(a) = sin (arctan(i;,)), iz (%o) is the ratio of top to bottom of the slope, æ is the slope of the road, 
and g is the gravity. 
Curve drag is calculated by (5): 


Ww, =o (5) 
With R is radius of the curve and A is coefficient determined by experiment, ranging from 450-800. With a 
rail gauge of 1435 mm, A=700; with a rail gauge of 1000 mm, then A=800. 


3. DESIGN TRACTION MOTOR DRIVE CONTROL 

The control design of the traction motor by FOC method [21] with three control loops: the current 
loop Ri, the flux loop Ry, and the speed loop R,, demonstrated in Figure 2. In the FOC control structure, the 
train speed is always greater than the rated speed of the traction motor, it is necessary to control the speed 
within the weak-field region. 
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Figure 2. FOC for traction motor 
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3.1. Constraints of stator voltage and current in weakening field region 

The two current equations of (1) show nonlinearity through the product of state variables isd, isq 
with input variables @,: Wsisq, Wsisq- Nonlinear structure features can be completely solved by nonlinear 
control methods. Therefore, in the control structure of the traction motor drive system, the bacstepping 
method chosen to design the speed and flux loop, and 7g in low flux region will be performed to ensure the 
speed control higher than the rated speed. 

Because the operating speed of the train is often greater than the rated speed of the traction motor, 
speed control is considered in two regions of rated flux and weakening flux corresponding to two speed 
regions W < Wygteq and W > Wrateq With the constraints of modules of stator current, voltage [21]: 


| 2 2 
Usa + usq < Usmax 


(6) 
«2 W 
: lsa + tog < I smax 
Stator voltage in steady state: 
Fa = fisa — WsOLsisg (7) 
Usq = Tsisq + WsLsisa 


The operating states leading to overvoltage all occur in the region of high frequency, so voltage drop 
on the resistor in (7) is smaller than that on the inductance, so it can be ignored. Substitute (7) in (6) to get: 


Sa isq 
Tma t T z (8) 
t) (C) 


In (8), drawing limitations of current and voltage in weakening field region in Figure 3. In Figure 3, 
green line is current circle, red line is voltage circle, black dashed line is the torque of the motor. It can be 
noticed, the higher the speed, the smaller the current isa and the motor torque will also be reduced, this is 
shown by the x-intersections. Figure 4 also indicated limitations of voltage us and currents is. 


Pra_rated 


speed 


Isq 


Tsa_ratea 


= | a > speed 
Rated field @ ateWeakening field 
region region 
Figure 3. Limitations of current and voltage in weakening Figure 4. Rated and weakening field region 


field (Wrated < w1) 


Rated speed range response in rated flux region: in this speed range even when the current is 
limited, the stator voltage never reaches the limitation. Under defined stator current module conditions, the 
two currents isd and isq must be controlled so that the working point is always at the maximum point of the 
torque- slip frequency curve. 

Torque equation: 


3 Li 


My = >Z —— 
M 2°? Letlrg 


İsa isq (9) 
From (9) my is maximum when isa=is4 and W, = Wrrateq = const 
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Lav = =—— (10) 


V3Ws5nLs(Isan) 


Speed range response in weakening flux region: in the flux attenuation range, there is an additional 
constraint of module of the stator voltage vector and the torque equation is calculated: 


m = 25 Lin Unax wrTr (11) 
M — 2°P Lye R2 (wrTrtwsTs)?+(1-ws0TswrTr)? 
. @y V3U „max LmN 
Pia = 0,9io amar am (12) 
(a) Un Lm 


The flux characteristic in (12) is inversely proportional to the rotor angular speed, which has adapted 
to the fluctuations of the Udc-link voltage and to the main inductance. 


3.2. Design flux and speed controller 

The backstepping technique is a recursive design method applied to the system which is comprised 
of some nonlinear subloops. The control law and Lyapunov control function are based on the state equation 
in each control loop to generate the virtual control signal for the next loops. From (1), using the backstepping 
method synthesizes usa and us, voltage signals to the PWM. 
The set values of speed and flux are w; and w;,, so the errors are defined as (13): 


Cy = Wr — Wr 
A 13 
Sa a 
Derivative of the above system of: 
Èw = O7 — Wr 
jeno (14) 
P = Pra ~ Pra 
Choose Lyapunov function V7: 
2 2 
V, = catep (15) 


2 


Derivative V7: 
Vy = eo (Wp — àr) + ey hra — ra) (16) 
For the speeds and fluxes to approach the set value, the derivative of V; must be negative and take the form: 
V, = -ke2 — kze% (17) 
Then the virtual control signals for the current are determined as (18): 
m Topai 1; 
lsa = 2Lm Wra + z Pra + kz2€ọ) 


2JLr 
32z2LmWrd 


(18) 


isq = (w; +m, + kieu) 


Where k; and kz are positive coefficients. 
The virtual control signal value is the set value for the next control loop. The errors of the currents 
are written as (19): 


Cin ~ isa zg isa (19) 
Cig = bsq T isq 
The derivative of the system of (8) combined with (1): 


r an r s% 1 1-o,. P 1-o 1 1-0 1 1 
êi, = lsa — lsa = lsa — Gr tsa t wslsqa t- — =o tu 

iq sd sd sd a m sd S°Sq oT, Lm Pra a Lin Pra Ls sd 20 
1 + i öz + yp EES 1 Y + 1 ii ( ) 
oT, Sq s*sd oT, Tq Lm rd Ls Sq 


soo L; LL 
Cig _ lsq lsq y lsq r a 
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Combining (7) and (8): 


Ted y 4% 
Cig = ae Wra + Yra + key) — lsa 
2Lm Tr 
21) 
a ely si Ze ; ( 
ei = ini OP +m, + kieo) — isq 


Choose the Lyapunov function V2 as (22): 


2 2 2 2 
esteyte; teig 


V = ————_. (22) 


2 


Taking the derivative of V2: 


x 1. 
= —k;eġ — key, = kzef, = kef, + eig (k3eig + tsa — (Gas arba t +@ al in Pra + 
-0 1 1 ok 1-o 
a in Pra P GLs Usa)) + Cig (k4ei; + lsq — (Cas =o) sq ® sisa + OT, A a Tw Ta Fi ne Ls l usq) 
(23) 


Where k; and k4 are positive coefficients. 
For the Lyapunov function to have a negative value, from (12): 
$ 1 1-0, f 
k3eig + isa — Gz + ae + ee ra + rq + JI “sa)) = =0 Kye, + isq — 
1 1-0, ; 1-0 
(c + zr isa — Wslsq + ar Vra age = Se alsa) = = 0 (24) 


oTs 


From (13) determine the control voltage signals of the system: 


. 1 1-o,. . 1-o 1 t-g 1 
Usa = OL (k3e;, + isa a (ele grisa + Wslsq + OT; Lm Dra + 5 Uyo Wrq)) Usg = 


ey 1 1- "1 1- 1 
OL (Ky4e;, + lsq — (Ce + -Diq w sisa + ae Pra B Zw Lin Wra)) (25) 


OTs oTr 


Where k; and k4 need to be corrected to ensure that the kinetic characteristics of the current converge faster 
than those of the speed and the flux. 


4. SIMULATION RESULTS 

Simulation parameters of train, traction motors which are collected from Metro line Nhon-Hanoi 
Station, Vietnam including 12 stations with 12.5 km in length in which 8.5 km elevated and 4 km 
underground have been shown in Tables 1 and 2. 


Table 1. Parameters Table 2. Parameters of traction motor 

Parameters of train Unit Quality Parameters Unit Value 
Mass of train at full load (M) kg 192000 Rated power (P) kW 185 
Numbers of motors (N) 12 Rated voltage (Un) V 525 
Maximum speed (vmax) km/h 80 Number of pole pairs 2 
Base speed (vb) km/h 40 (Zp) 
Maximum acceleration m/s? 1 Rated speed (n) revolution/minute 1485 
Coefficient a 0.0115070 Inertia moment (J) kg/m? 4.2 
Coefficient b 0.0003494 Stator resistance (Rr) N 0.0127 
Coefficient c 0.00005497 Rotor resistance (Rs) N 0.0127 
Wheel diameter (Dwh) m 0.84 Stator inductance (Ls) N 0.0109 
Transmission ratio (i) 5.67:1 Rotor inductance (Lr) a 0.0109 
Gearbox performance (Nmech) 0.85 DC-link voltage V 1000 
Motor performance (Nem) 0.95 (Udc) 
Inertia moment of train (Jeq) kg.m? 87.67 Total inertia torque (J) kg.m? 91.87 


Simulation results were conducted for trains running between two stations: Kim Ma-Cat Linh with 
the distance of 1508 m including 3 phases: accelerating, holding, and accelerating. In which, running time of 
92.2 s, train speed up to 45 km/h: acceleration time t=15.46 s, distance S=99.84 m; holding running time 
t=65.68 s, distance S=1337.46 m; braking time t=11.06 s, distance S=70.7 m in Figure 5. 
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Speed curve from Kim Ma to Cat Linh station 
| | | | — 


v(km/h) 
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Figure 5. Response of speed curve 


Figure 6 shows response of consumption power: the train mobilizes the maximum capacity from 0 
to 13.22 s; from second 13.22" to second 15.46" having the constant power, hyperbolic torque corresponding 
to the motor operating in the weakening flux region; from second 15" to second 81", the power consumption 
decreases corresponding to the train running at steady speed. After second 81", the negative power shows 
that the energy is returned to the grid corresponding to the train running during the braking. 
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Figure 6. Consumption power of a traction motor 


Figures 7 and 8, from second 13.22 th to second 15.46 th shows the period of weakening flux and 
torque when motor speed is greater than rated speed. This is a distinctive feature of electric train loads 
compared to other types of loads in the industry. In the holding phase, the torque decreases, while braking 
phase, negative torque indicates the regenerative braking process, where energy is returned to the grid. 
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Figure 7. Torque of a traction motor 
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Figure 8. Response of flux 


5. CONCLUSION 

When urban trains operate during the acceleration phase at speeds greater than 40 km/h, it is often 
preferred to control the traction motor in the weak-field region to increase the train's speed while still 
ensuring maximum torque control. Therefore, this paper proposes to control the speed of traction motor in 
weakening field region with the constraints of stator voltage, current modules, using backstepping technique 
with simulation parameters collected from Metro Nhon—Ha Noi Station. The simulation results have shown 
the responses of train speed, flux, power, and torque are completely consistent with 3 train operating modes: 
accelerating mode, holding mode, and braking mode. 
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